Generation of Spin Entanglement in Nonequilibrium Quantum Dots 
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We propose schemes for generating spatially-separated spin entanglement in systems of two quan- 
tum dots with onsite Coulomb repulsion weakly coupled to a joint electron reservoir. An enhanced 
probability for the formation of spin entanglement is found in nonequilibrium situations with one 
extra electron on each dot, either in the transient state after rapid changes of the gate voltage, or 
in the steady state with applied bias voltage. In both cases so-called Werner states with with spin 
singlet fidelity exceeding 1/2 are generated, which indicates entanglement. 

PACS numbers: 03.67.Mn,73.23.Hk,73.21.La,73.63.Kv 



I. INTRODUCTION 

The entanglement of quantum states is one of the cor- 
nerstones of quantum information processing^ Entan- 
gled photons have been used in experiments in quan- 
tum communication and cryptography.— For electrons in 
a solid-state environment recent progress has been linked 
to advances in fabrication technology for nano-scale de- 
vices^^ The availability of an electron spin entangler in 
a solid-state environment would allow the implementa- 
tion of quantum communication schemes with electron 
spins££ Several schemes have been suggested for the pro- 
duction of spatially-separated entangled electrons in solid 
state systems. Many of them rely on extracting the en- 
tangled electrons of a Cooper pair from a superconductor 
and separate them into two normal leads, 7 Luttinger liq- 
uids,— or to two leads through two quantum dots.— Others 
are based on separating the electrons forming a spin sin- 
glet on a double-quantum dot,™ using interference effects 
in a quantum dot in the cotunneling regime,— separat- 
ing a pair of entangled electrons from a singlet state by 
a triple quantum dot ^2. or scattering off magnetic impu- 
rities M- 

In this article we show that a pair of entangled elec- 
trons can be created by driving out of equilibrium a sys- 
tem of two quantum dots with onsite Coulomb repulsion 
and weak coupling to a joint electron reservoir. Specifi- 
cally we consider the two setups depicted in Fig.[TJ Elec- 
trons enter the dots from the reservoirs, and we consider 
the nonequilibrium state with one electron on each dot. 
In setup a) we study the transient behavior after quickly 
pushing the dot levels from above to below the Fermi 
energy of the lead, and find an enhanced probability for 
the singlet state as compared to a triplet. In setup b) 
we drive the system out of equilibrium by applying a 
bias voltage between left and right leads. Depending on 
the polarity of the applied bias, we find in the steady 
state an enhanced probability of either the singlet or the 
triplet states. The mixed states with two electrons in 



FIG. 1: The setups: a) Two quantum dots (u and d) are 
coupled to a joint electron reservoir (L). b) In addition to a), 
the quantum dots are coupled to two independent reservoirs 
(Ru and Rd) on the right. 



the two dots represent so-called Werner states*^ In the 
case where the electrons entered from the common (left) 
reservoir we find regimes where the Werner fidelity is 
larger than 1/2, which implies a high probability for the 
formation of a singlet state. 



II. MODEL 



The Hamiltonian of the system is 

" = "dots + -ff leads + "tunnel- 



(1) 



The two quantum dots, i = u, d (up and down), described 
by 



-Hdots — 



Si C, 



(2) 



contain each a single, spin-degenerate energy level Si. In 
general the dot levels are detuned by Ae = e u — Ed- We 
assume strong Coulomb repulsion within each dot, U ^> 
&bT, eV, r which suppresses double occupancy of each 
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dot. (Our analysis can be generalized to finite interdot 
charging energy, which does not change the conclusions 
qualitatively.) The leads 

-Hleads = ^ ^ £rk a tka a rka , (3) 
r ka 

with r = L,R u ,Rd, serve as equilibrium reservoirs with 
electrochemical potentials fi r . The tunneling between 
leads and dots is modeled by 

-fftunnol = ^2 ^ (tri c\ a a rk<y + h.C.J . (4) 
r kai 

The tunneling strength of quantum dot i to reservoir r is 
parametrized by T ri = 1-Kt 2 ri N r , where N r denotes the 
reservoir density of states, and we chose a gauge in which 
all the tunnel amplitudes t ri are real. We further define 
r. r = J2i F-ri/2. For setup a) we have Tr u = T^ d = 0. 

Our proposal is based on the observation that the 
states in the common left lead are only coupled to a cer- 
tain linear combination of the two quantum dot states. 
If Coulomb interaction were absent, U = 0, filling the 
double dot with two electrons with opposite spin from 
the common reservoir would lead to the product state 

(*Lu4r + *Ld4r) (*Lu<4 + *Ld4») 1°) - 0) + 

^Lu^Ld (|f, o) — |<7, o)) + i| d |0, oo). For strong Coulomb 
repulsion, however, the parts that involve double occu- 
pancy of either dot are projected out, and the final state 
is, \o, o) — \o, o), no product state but a spin singlet. No 
triplet component, although energetically degenerate to 
the singlet, is generated. 

In realistic situations various mechanisms will relax 
the imbalance between the population of spin singlet and 
triplet states, e.g., tunnel coupling to the right reservoirs 
shown in Fig. [TJa), or a finite detuning Ae. Furthermore, 
a coupling to an external bath which mediates spin-flip 
processes or creates a phase difference between the dot 
states causes an equilibration between singlet and triplet. 
In this article, we study in detail nonequilibrium scenar- 
ios characterized by the competition between the creation 
of singlet and triplet states and the relaxation. 

III. KINETIC EQUATIONS 

For this purpose we employ the real-time diagrammatic 
technique developed for single quantum dots^ and ex- 
tended to multi-dot systems! 16 ' 17 In this technique the 
electronic degrees of freedom of the leads are integrated 
out, which results in an effective description in terms of 
the degrees of freedom of the dot subsystem only. The 
dynamics of the latter is then described by a reduced den- 
sity matrix with elements = (\x)(x'\/> where \ an d 
X 1 label the double-dot states, and (. . .) denotes quan- 
tum statistical expectation values. In the present case, 



the Hilbert space of the quantum-dot degrees of freedom 
is spanned by 9 basis states |x u >Xd), with Xi € {AIM} 
denoting the occupation of dot i. 

The time evolution of the reduced density matrix in the 
Markovian limit is governed by the kinetic equations 15 

J t vt + i{E x . - E x )p£ = W tf $ ■ ( 5 ) 

x"x"' 

The energy difference E x > — E x between states x' an( i X 
leads to a time-dependent phase of the off-diagonal ma- 
trix elements. Transitions due to the tunnel coupling to 
the leads are described by the kernels W£Jt> , the general 
form of which are given in Refs. [Hill. In the following 
we restrict our attention to the limit of weak coupling 
and small detuning Ae, where it is sufficient to evaluate 
the kernels in first order in the tunneling strength T r i 
and to zeroth order in Ae. 

To proceed it is convenient to switch to a basis {|x)} 
which reflects the symmetries of the problem. One of the 
basis states is the empty-dots state |0) = |0, 0). For two 
electrons, one in each dot, the natural basis states are 
the spin singlet |S) = (| |,|) - | |,f))/V2 and triplet 
states |T+> = | T,T), |T ) = (| U) + | U»/\/2, and 
|T_) = | 1,1). The states with one electron in the dou- 
ble dot can be characterized by the physical spin o of 
the electron, as well as by an isospin defined in the 2- 
dimensional Hilbert space spanned by the two orbital dot 
levels. One natural quantization axis n for the isospin op- 
erator I CT is the one in which the eigenstates of \ a ■ n are 
|+)i CT . n = W, 0) and |— )i o .. n = \0,o), corresponding to the 
electron in dot u and d, respectively. This is motivated 
by the observation that both the Coulomb interaction 
and the detuning Ae in the Hamiltonian is diagonal in 
this isospin basis. An alternative choice is the axis m 
defined by |+)i <r . m = (t Lu \o, 0) + t Ld |0, a)) / ^ij n + ij d 

and |-) lCT . m = (i L d|M)-tLuM)A/C+ t 2 hd . This 
reflects the fact that filling the double dot by tunnel- 
ing with one electron from the left lead generates the 
isospin component |+)i„. m onlyJ£ In this sense, the left 
lead can, in analogy to magnetism, be viewed as a fully 
isospin-polarized lead with only + isospin-electron states 
available. The right reservoirs in setup b), on the other 
hand, correspond to an isospin-unpolarized lead. In gen- 
eral, the two axes n and m are not orthogonal, except 
for the symmetric case when the tunneling strengths to 
dot u and d are equal, as can be seen from n • m = 

(r L u-r Ld )/(r Lu + r Ld ). 

The total Hamiltonian is invariant under rotations in 
spin space, i.e., spin is a conserved quantum number. 
Spin symmetry implies (If) = (Ij) = 1/2 as well as 
Pt_ = Pt = Pt + = Pr/3, which reduces the number 
of independent density matrix elements. Those are the 
isospin I and p = (po,Pi,Ps,Pr), where pi = Yl.iaVia is 
the probability for single occupation. In this representa- 
tion, the kinetic equations read 



3 



/-4/ r l-/r 



P = 



r=L,R 



dt 



= r r 



4/,- 


-1 - fr 





fr/2 


V o 


3/r/2 


»+ (h 









2 - 2/ r 
-2 + 2/ r 






2 - 2f r 


-2 + 2/ r 



P + Ti 




(I-m) + 2r L / L _ x (I-n)(m-n) 



m + r L 



^pi - 2(1 - / L )p s 




r r (l + / r )I + Ae(nxI), 



(6) 



r=L,R 



where / r = [1 + cxp(/3(e — /v))P is the Fermi dis- 
tribution of the electrons in lead r. Here we intro- 
duced, apart from the detuning Ae also the average 
e = (e„ + £d)/2 of the dot energies. The level de- 
tuning is renormalized by the tunneling and given by 



Ae = Ae - 



r L „ 



In(f2)-B** 



2 + ^ 2tt 



Ml) 



where 13 is an high-energy cutoff provided by either 
Coulomb interaction U or bandwidth of the leads. 



IV. RESULTS 

A. Spin Entanglement in Transient States 

Inspection of Eqs. ^ reveals how an imbalance of sin- 
glet and triplet states can occur. When filling the empty 
double dot with one electron, a finite isospin along m is 
generated. This in turn, blocks the generation of triplet 
states as opposed to singlet states when filling the dou- 
ble dot with a second electron. This mechanism becomes 
most transparent for Tl u = lYd, Fr u = ^R d = and 
Ae = 0. In this case, the two equations 

j t PT = 3r L / L (a _ i . m ) _ 2 r L (l - f L ) PT (7) 



d_ 

dt 



-3r L / L (y - I - mj+2r L (l-/ L )p T , 

(8) 

decouple from the rest. This motivates proposal a) for 
generating spatially-separated spin entanglement. If we 
prepare the system in an empty state (by tuning the dot 
levels well above the Fermi energy of the lead) and, then, 
suddenly push the dot levels well below the Fermi energy 
of the left lead, — e ^> keT, Tl, the double dot will be 
charged with two electrons that form a spin singlet, while 
no triplet component appears. The time dependence of 
the singlet generation is illustrated in Fig. [2] 

Coupling to an external bath, which flips the spin of an 
electron or generates a relative phase between the | f, [) 
and | |,t) states, induces relaxation from the singlet to 
the triplet state. To model these processes we introduce 
phenomenologic relaxation rates Tg— >t ) rs^T ± , Ft ^Sj 
Ft ± ^s> rTo^T ± , and Ft ± ^t - To be specific, we choose 
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FIG. 2: Upper panels: time evolution of the probabilities for 
a singlet and a triplet state. Lower panel: the corresponding 
Werner fidelity. For the perfectly symmetric setup, Ae = 0, 
Flu = Ld, in the absence of spin relaxation, Fs^t = 0, 
curve i), we find F = 1. The Werner fidelity is reduced for 
either ii) nondegenerate dot energy levels, Ae = Fl, iii) asym- 
metric coupling TlcI = 0.1 Tm, iv) a finite spin relaxation rate 
Fs^t = 0.2 Fl- The high-energy cutoff is set to D = 100 k B T. 



all of them to be equal, such that we get an effective tran- 
sition rate Ts^t = which conserves the symme- 
try between the triplets, pt ± = Pt - (A different choice 
of these parameters does not change the conclusions qual- 
itatively.) Furthermore, a finite detuning Ae and/or fi- 
nite asymmetry of the tunnel couplings Tl u TLd, lead 
to a mixture of singlet and triplet states, producing a 
Werner state« described by the density matrix 



W{F) =F|S)(S| + (1-F) 



|S)(S| 



(9) 



The parameter F defines the Werner fidelity. Werner 
states play a crucial role in entanglement purification 
protocols ) 19 ' 20 and the Werner fidelity gives a convenient 
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FIG. 3: Time evolution of the total charge during the dis- 
charging of an initial singlet state compared to an initial 
triplet state. The system is assumed to be perfectly sym- 
metric, Tm = TLd, Ae = 0. Starting from a singlet the 
system empties quickly, but it remains singly occupied if we 
start from a triplet state. 
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FIG. 4: The stationary Werner fidelity F vs. bias voltage 
for e = and different ratios of the coupling strengths, 
Tr/Fl = 2, 1,0.5. The bias is applied symmetrically, — fit, = 
eV/2 — /Ur u = /iR d . The inset shows the corresponding sta- 
tionary overall probabilities ps +Pt to find the system doubly 
occupied. 



measure for the possibility to extract entangled states 
from a set of Werner states by such protocols. It has 
been further shown that for Werner fidelity 1/2 < F < 1 
there exist purification protocols to extract states with 
arbitrary large entanglement whereas for F < 1/2 the 
Werner state has to be considered as unentangled. 

Solving the kinetic equations for the reduced density 
matrix for system a) we see that Werner states with fi- 
delity F = ps/(ps + Pt) > 1/2 are accessible also for 
asymmetric tunneling, detuning and finite spin-flip relax- 
ation, see Fig. [2] For weak detuning Ae the probability 
to generate a triplet scales with pt ~ (A£/2Fl) 2 . 

To create and detect an enhanced spin-singlet fidelity 
and to measure the relaxation time between singlet and 
triplet we propose the following scheme that is similar 
to the experiment performed in Rcf. |3|. (i) Prepare the 
system in an empty state, (ii) Push quickly (i.e. on a 
time scale faster than both the relaxation times for the 
isospin-polarized state and for the singlet-triplet transi- 
tions) the dot levels down well below the Fermi level. 
As explained above the double dot will preferably fill up 
with two electrons forming a spin singlet, (iii) Wait some 
given time T. As a function of T, the imbalance between 
singlet and triplet decays exponentially on the time scale 
given by the relaxation rate, and the Werner fidelity is 
reduced, see Fig. [5J 

To prove that the obtained state, indeed, has an 
enhanced Werner fidelity, we analyze how the double 
dot is depleted. Depending on whether the initial 
state is a singlet or triplet, it is possible or impossi- 
ble to extract the two electrons by tunneling to the 
common left lead. This can be seen by realizing that 

(^LuCuo- + t Ld Cdcr) (thuCua + *LdCdff) (|cr, d) T I*, a)) = 

(1 ± l)^Lu^Ld|0) is finite for the singlet but vanishes for 
the triplet state, i.e., only one of the two electrons form- 



ing the triplet can leave. As a consequence, the proposed 
protocol continues in the following way. (iv) Push the 
dot levels up well above the Fermi level quickly (again 
faster than the relaxation rate for the isospin-polarized 
state), (v) Wait some time larger than I/Fl but shorter 
than the relaxation time of the isospin-polarized state, 
(vi) Measure the total charge on the double dot. If the 
charge is zero then the doubly-occupied state was a 
spin singlet, whereas if the measured charge is one, it 
was a triplet. To illustrate this we show in Fig. [3] the 
total double-dot charge as a function of time for the 
two cases that the double dot initially accommodated a 
singlet or a triplet, respectively. The measurement of the 
total charge on the double dot could be performed by a 
close- by quantum-point contact. This does not introduce 
an additional relaxation mechanism for either the isospin 
or the singlet and triplet states as the quantum-point 
contact is only sensitive to the total charge. 



B. Spin Entanglement for Finite Bias Voltage 

Spatially-separated spin entanglement is found also in 
a steady-state situation in the setup b) of Fig. [TJ Here 
we consider the system to be driven out of equilibrium 
by a bias voltage between the left and the right side. 
To keep the discussion transparent we assume in the 
following symmetric couplings, Tl u = TLd = IY and 
rR u u = rR d d = Tr, equal electrochemical potentials in 
the right leads /ir u = /iR d , and vanishing detuning of 
the dot levels Ae = 0. The leads on the right hand side 
couple to all isospin components in the same way. In a 
magnetic analogue such a situation corresponds to a dot 
coupled to one ferromagnetic and one nonmagnetic lead 
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for which, at large bias voltage, spin accumulation occurs. 
Similarly, in the present model a finite isospin is accu- 
mulated in the double dot in the stationary limit. This, 
again, leads to an imbalance of singlet and triplet state 
probabilities. The polarity of the bias voltage determines 
whether the Werner fidelity is larger or smaller than 1/4. 
If the bias voltage is applied such that the double dot 
is charged from the left and decharged to the right lead 
the isospin polarization is in + direction, and singlets are 
preferred. In this regime the Werner fidelity saturates at 
F = (3 T L + 2 r R )/(6 T L + 2 r R ) which goes from 1/2 for 
Tl S> Tr to 1 for Tl <C Tr. We have to remark that the 
fidelity approaches 1 only linearly for IY <C Tr, whereas 
the overall probability to find the double-dot system dou- 
bly occupied vanishes quadratically ps +pt ~ 2 (IY/Fr) 2 
at the same time. If the bias voltage is applied in the op- 
posite direction, triplets are more likely. 

V. CONCLUSIONS 

For an experimental realization of our proposal one 
needs to coherently couple two quantum dots to a joint 
reservoir, as has been demonstrated e.g. in Ref. [U The 
spatial separation of the two dots is only limited by the 
phase-coherence length, which can be several micrometer 
in typical semiconductor structures. The formation of 
an enhanced spin-singlet fidelity requires tunneling rates 
larger than the spin decoherence time. Reported val- 
ues^ 2 - of T 2 * of the order of 10 ns correspond to a lower 
limit of r of the order of /ieV. For tunnel couplings T 
larger than k^T higher-order processes such as cotun- 
neling and Kondo-assisted tunneling become important. 



These are neglected in our quantitative analysis but they 
do not change our prediction qualitatively. In fact, for 
symmetric tunnel couplings the Hamiltonian acquires a 
block structure and the Hilbert subspace containing the 
triplet states decouples completely from the one for the 
empty double dot. In conclusion, the experimental real- 
ization of our proposal should be feasible by nowadays 
technology. 

In summary, we proposed two schemes for the gener- 
ation of spin entanglement between two spatially sep- 
arated electrons in a double-dot system driven out of 
equilibrium. The underlying mechanism is fundamen- 
tally different from those that rely on a singlet-triplet 
energy splitting, where entanglement is generated by a 
relaxation of the system to the spin-entangled ground 
state. In contrast, we suggest two schemes in which en- 
tanglement is a consequence of a coherent coupling of two 
quantum dots to one common lead in combination with a 
strong onsite Coulomb interaction to prevent double oc- 
cupancy of each individual dot. We emphasize that our 
proposal does not require a finite singlet-triplet splitting. 
The quick formation of the entangled state on a time 
scale given by the tunneling instead of a singlet-triplet 
relaxation rate, may be an advantage in the context of 
quantum information processing. 
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SFB 491, the Swiss SNF and NCCR Nanoscience, as well 
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